Detecting drift of quantum sources: not the de Finetti theorem 
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We propose and analyze a method to detect and characterize the drift of a nonstationary quantum 
source. It generalizes a standard measurement for detecting phase diffusion of laser fields to quantum 
systems of arbitrary Hilbert space dimension, qubits in particular. We distinguish diffusive and 
systematic drifts, and examine how quickly one can determine that a source is drifting. We show 
that for single-photon wavepackets our measurement is implemented by the Hong-Ou-Mandel effect. 



Ever since its first experimental implementation in 
1993 PJ, quantum tomography (or, more generally, 
quantum-state estimation [2,) has become an important 
tool in the field of quantum information science (for re- 
views, see [Sill])- From the results of different measure- 
ments on many instances of a quantum system, one infers 
a density operator, po (plus "error bars") that describes 
best (in some more or less well-defined sense) the state 
of each instance. One may well wonder why one assigns 
just one density operator. A crucial role in this context 
is played by the de Finetti theorem (for details, on both 
the infinite and finite versions of the theorem, see [3 [6]): 
if one has an extendible permutation-invariant sequence 
of N quantum systems, then one can assign a quantum 
state of the form 

= J dpP{p)p^^ (1) 

to the collection of N systems, with P[p) a probability 
density over density operators. Quantum tomography 
then succeeds in making the distribution P{p) more and 
more narrow, sharply peaked around some po- In fact, in 
the limit of TV — >• oo, one has P{p) — 5{p — po). 

We are interested here in the case where the assump- 
tion of permutation invariance does not hold, and where, 
consequently, the de Finetti theorem does not apply. The 
most relevant case is that of a (slowly) drifting source. 
For example, it is well known that a laser displays phase 
diffusion: when one considers two light pulses emitted by 
the same laser with a short time delay r between them, 
there will be a (random) phase difference whose aver- 
age magnitude increases with r. Of course, even in this 
case, one could average over all emitted light pulses, say 
TV instances, to arrive at a single average density ma- 
trix. Indeed, if done correctly the averaging procedure 
restores the permutation invariance, but (i) the average 
density matrix depends on the number N, and (ii) the 
averaging procedure throws away potentially useful infor- 
mation. For example, if we are interested in the purity 
of our quantum states, the single state estimate will be 
too conservative. 

In this paper we set ourselves the task of figuring out 
how one could detect whether (and how) a source is drift- 
ing. In principle, for detecting drift one could still use a 
variant of quantum tomography: for example, we split 



our N quantum systems into two groups of size N/2 
each: the first half (chronologically) and the second half. 
For each we estimate a single density matrix: and if 
the difference between the two estimates is (not) statis- 
tically significant then we conclude our source is (not) 
likely drifting. This method works to some extent, but is 
still subject to the same two objections mentioned above. 
Moreover, it has been known for a few decades that for 
the detection of given physical quantities (such as a par- 
ticular matrix element of the density matrix) a targeted 
method is always superior to performing full tomogra- 
phy Therefore, we propose and analyze a different 
method directly targeted at detecting drift. A difference 
with the above-mentioned method [7| is that we consider 
a quantity determined by pairs of density matrices. 

Consider what one would measure to detect phase dif- 
fusion of a (pulsed) laser in the special (but relevant) case 
where one assumes the laser pulses can be described by 
coherent states with some fixed (and known) amplitude 
but a diffusing phase (relative to some phase standard). 
One would take pairs of the output laser pulses, and split 
them on a 50/50 beamsplitter in such a way that one par- 
ticular output would be the vacuum if their phase differ- 
ence, S(j), would be zero. That output's intensity is then 
/ = — exp(z(5(/))p/2, if \a\ is the amplitude of the 

laser pulses. Thus measuring this intensity determines 
the phase difference directly. 

Now how do we generalize this measurement to arbi- 
trary quantum systems (in particular, qubits)? We first 
note that the intensity / can also be written in terms 
of the overlap between the two input states, call them 
p = \a){a\ and p' — \a'}{a'\, since exp(— 2J) = Tr{pp') 
in this case. Thus, our choice of generalization will be to 
measure the overlap between pairs of instances of quan- 
tum systems from one and the same source. In other 
words, we propose to measure the swap operator V, de- 
fined in terms of basis vectors {\i)} and {\j)} of the two 
(isomorphic) Hilbert spaces of two instances numbered 
m and n from our source by 

V^J2J2\')r,.{j\®\j)rM (2) 
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The expectation value of V equals the overlap 

Tr(p„ (g) PnV) = TT{pmPn)- (3) 
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(Note the left-hand side contains the tensor product of 
two density operators, the right hand-side their matrix 
product.) Each measurement of V yields one of its two 
eigenvalues, ±1, and so only after multiple measurements 
will one obtain a statistical estimate of the overlap. [And, 
as a bonus, if the two density matrices are identical, 
then this measurement in fact measures the purity [5J [5] , 
P = Tr(/9^).] By comparing the overlap between ad- 
jacent copies, where \m — n\ = 1 with the overlap be- 
tween outputs that are farther apart, |m — n| > 1, we 
obtain information about whether the source is drifting: 
if the source is not drifting, the overlap is independent of 
\m — n\. 

In order to infer more detailed information about the 
character of the drift or diffusion (beyond the mere 
statement that the source is or is not stationary), we 
need some simplifying assumptions about the sequence of 
states (the space of all possible output states of N copies 
is too large to be either measurable or tractable). Here 
we make the following two assumptions (one of which has 
been implicitly used already in the above description): 
(a) the states are independent, (b) the drifting process is 
Markovian, such that the overlap between two copies m 
and n only depends on |m — n|. So we write the state of 
N systems produced by our quantum source as a tensor 
product [12] 



P^^-' ^ PN® PN- 



Pl- 



(4) 



In this case, we get 



Vnra = Tr[p„ ® p^V] = \ [P„ + - ^Tr[A^„] , (5) 
where Pfe = Tr[(pA;)^] is the purity of system A:, and A„i„ 



is the difference between the two states m and n: 



A, 



(6) 



As special cases of nonstationary sources we consider 
both diffusion and systematic drift. Consider a (Marko- 
vian) process where 



UrPnUl- 



(7) 



Diffusive drift occurs when \Jr is a random unitary ma- 
trix, picked from some distribution; a systematic drift 
occurs when is fixed. In either case, the purity of 
Pn is independent of n: P„ = = Pi. We can take 
stochastic averages over the random distribution of uni- 
taries, which we will indicate by a bar, to get 



1 



= Pi--Tr[A2J. 



(8) 



In case the drift process is purely a systematic drift, each 
unitary Ur is the same, and we get 



Tr[A2J = |n-m|^Di 



(9) 



If the process that changes the states p„ is diffusive 
(for example, the random distribution of Ur is a Gaus- 
sian centered around the identity), then we get a linear 
relationship between the overlap and the distance |n— m|. 



Tr[A 
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(10) 



for some diffusion constant D2 |17) . In this case, mea- 
suring the swap operator between neighboring copies, for 
which |n — m| = 1 and on copies with |n — mj — 2 gives 
us both the purity Pi and the diffusion constant D2'. 



Pi = 2Tr[A2.„+J 



Tr[A2, 



D2 = Tr[A2^„+2] 



Tr[A2 



n,n+lJ 



(11) 

(12) 



(And similar relations hold when the drift is purely sys- 
tematic.) One way to check which sort of drift pro- 
cess one actually has, diffusive, systematic, or a com- 
bination thereof, is to measure in addition the quantity 
Tr[A^ n+sli ^-nd then calculate the ratio 



Tr[A2,„+2] 



- Tr[A2^„+i] 



Tr[A2,„+3] 



(13) 



If the ratio is 1, one has a purely diffusive process, if 
a = 3/5 one has a systematic drift, and in all cases in 
between one has both diffusive and systematic drifts. To 
see how the number a is determined when there is a com- 
bination of systematic and diffusive drifts, let us consider 
the simplest case of a qubit source. We model the drift 
process with a unitary matrix Ur = exp(i(5r'- a), with 
(5 ^ 1, (7 a vector containing the three Pauli matrices, 
and a random vector r that consists of both a diffusive 
part and a systematic part, 



r^pr const + (1 -p)7\ii 



ffusive 



(14) 



with a normally distributed random vector rdiffusivc and 
a constant (unit) vector rconst (and < p < 1). In this 
case a depends on p and on the ratio of the two constants 
Di and D2, with D i = (5^Tr[fco„st • S,pY' and D2 = 
(5^Tr[fdiffusivc ■ S,pY: 



p2 + ^iL(_5p2 ^ lOp- 5)' 



(15) 



The next question we consider is, given a source of quan- 
tum states, how quickly can we determine (by measuring 
the swap operator) whether the source is drifting? Let us 
first consider the case of pure diffusive drift. We could, 
for example, measure the swap operator between states 
that are 1 step and 2 steps apart, respectively, and see if 
the two numbers are equal or not. Suppose we write the 
overlaps as 



Tr(p„p„+i) =Pi -i?2/2, 
Tr(p„p„+2) = Pi - D2. 



(16) 



for some drift constant Di . 
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Suppose we have N data sets of the measurements of 
both Tr(p„p„+i) and Tr(y0„p„+2). We get the measured 
frequencies and , respectively, of the measurement 
outcomes ±1 in the two cases. The average values are 



Vi := Tr(p„p„+i) = /f - , 

V2 Tr(p„p„+2) = f+ -f^, (17) 

and the standard error bars are (for large enough N) 



AVi 



AV2 



ftfi 



N 



N 



(18) 



To decide that the source is drifting, the values of Vi 
and V2 should not overlap within their error bars. The 
necessary condition for that is 

^AFi + ^A^a <Vi-V2= D2/2. (19) 



frequencies in terms of Di and Pi 



From assumption (16) and Eqns. (171, we can write the 

D2. 



ft - + ^1 - /r = ^(1 - ^1 + 
ft = \{i + Pi- D2); f^ = i(i - Pi + D2) 



(20) 



and using all this in Eq. ( 19 ) we can solve for the mini 



mum necessary number of measurements 



^mi„ - [-^(^ (1 + Pi - ^)(1 -Pl + ^) 

+ ^{l + Pi-D2){l-Pi+D2))f. (21) 
To detect a systematic drift, a similar calculation gives 



A^.. = [3^(V(l + A-^)(l-Pi + ^) 

+ v/(l + Pi - 2Z?i)(l - Pi + 2i?i))]2. (22) 

The number of measurements needed, for both diffusive 
and systematic drifts, is depicted in Fig. [T] for various 
values of Pi . 

In principle, one could detect a drifting source a lot 
faster if one measured the swap operator on states that 
are k > 2 steps apart (in addition to measuring states 
1 step apart), simply because \Vk — Vi\ will be larger. 
In an actual experiment, however, the larger the dis- 
tance between two copies, the longer the earlier copy 
would have to be stored in memory. We could model 
the decoherence that the earlier copy undergoes as fol- 
lows: assume that there is a typical decoherence time 
scale r, which, e.g., drives any state towards the totally 
mixed state. That is, if we keep a system for time i, 
then (exp(-t/T)p+(l-exp(-t/T)l/D, with £» the 
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FIG. 1: How many measurements do we need to figure out 
that a source is drifting? Obviously, the larger the drift is 
(as measured by the parameter D2 for diffusive drift or Di 
for systematic drift), the fewer measurements we need. Top; 
systematic drift, bottom: diffusive drift. 



dimension of the Hilbert space of our quantum system. 
Then assume that the time needed to produce one copy 
is er with some (hopefully small) number e. Then we can 
write the overlap between states n and n -I- fc as 



Pfe = Tr(p„p„+fe) 

= e-*'-^Tr(p„p„+fe) + -^(1 - e-'=^)Tr(lp„+,). (23) 

The inferred overlap between copies n and n + k follows 
from the measured Pfc by multiplying it with exp(A:e) 
(and subtracting a known quantity): so the error bar in 
the overlap multiplies by the same number. This error 
thus becomes substantial once ke becomes of order unity, 
so that is where we would expect the method to use copies 
a distance k apart to break down. In Fig[2]we plot the 
number of measurements needed for various values of e 
as a function of k, and we can indeed see that for too 
large values of fc, the required number of measurements 
increases exponentially with k. The optimal k, of course, 
depends on the specific decoherence process, and also on 
how fast the source is drifting, but seems to be around 
fee « 4 in our example. 

Finally, we wish to note that in the case of two inde- 
pendent single photons, when they are viewed as quan- 
tum systems with an infinite-dimensional Hilbert space 
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FIG. 2: This plot shows the minimal number of measure- 
ments needed to detect drift, as a function of the (temporal) 
distance between states measured, k, for various values of the 
decoherence parameter e, for Pi = 1, D — 2 and D2 = 0.01 
(see main text for definitions). The minimum of each curve 
determines the optimal distance k between states to be mea- 
sured (in addition to distance- 1 overlaps). 

describing polarization, spectral, and transverse spatial 
degrees of freedom, the swap operator can in fact be mea- 
sured via the Hong-Ou-Mandel interference effect [TU] . 
This can be shown as follows. We consider two sin- 
gle photon wavepackets impinging on two different input 
ports (denoted A and B) of a 50/50 beamsplitter. We 
write the two (mixed) input states in terms of creation 
and annihilation operators and a (for port A) and 
and b (for port B) as: 

PA^j:kiPkial\0){0\ai (24) 

PB^En^1nmbi\0){0\b^, (25) 

where the subscripts stand for the mode properties (po- 
larization, frequency etc.) other than their propaga- 
tion direction. The combined input state is then pin = 
Pa®Pb- This state gets transformed by the 50/50 beam- 
splitter in the following way: 

Pout = Efe,„™ + *4)(*4 + 4) (26) 

|0) (0| {ci - idi){-icm + dm), 

where c and d now denote operators of the two output 
ports C and D. To get the probability Pec of getting a 



coincidence count, i.e., photo detections at both output 
ports C and D, we take a partial trace: 

Pec =E„(l.lc(l^ldPout|l.)clUd (27) 
^Ers {0\CrdsPout4dl\0) 

This simplifies to 

Pec ^ ^^Pkk^ qnn - ^ ^PklQlk- (28) 
k n kl 

The first two sums are the traces of the density matrices 
and therefore equal 1. It is easy to see that 

TripAPs) = ^ PkiqnmSnlSmk = ^PklQlk, (29) 

klnrn kl 

SO that we get the simple relation 

2Pec = l-Tr(pAPB). (30) 

Thus, as announced, the HOM effect measures the over- 
lap between two input states, and hence the swap oper- 
ator. (And so, if the two single-photon input states are 
identical, then the HOM interference measurement mea- 
sures the purity of the input states. Note that this is 
different from the measurement of single-photon (spec- 
tral) purity implemented recently in Ref. , which also 
makes use of the HOM effect, but with a known coherent- 
state input in the other input port.) Of course, the HOM 
effect has been measured many times in the context of 
characterizing single-photon sources (see, e.g., [HI [13]), 
but never, as far as we know, systematically on copies 
more than the minimum distance apart. We also note 
that for the polarization degree of single photons, the 
overlap has been measured [T3], following ideas from [H]. 

In conclusion, we proposed the measurement of the 
swap operator as a means to detect the drifting of a 
quantum source. This measurement complements quan- 
tum tomography, which produces an estimate of a single 
average density matrix, by partially characterizing how 
this estimate would change over time, for instance, dis- 
tinguishing between diffusive and systematic drifts. We 
also analyzed how many measurements are needed to 
determine that a source is drifting, including the influ- 
ence of decoherence on the precise measurement strategy. 
We showed the swap measurement on pairs of single- 
photon wavepackets is implemented simply by the Hong- 
Ou-Mandel effect. 
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